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Abstract. Measurement of a quantum system provides information 
concerning the state in which it was prepared. In this paper we show how 
the retrodictive formalism can be used to evaluate the probability 
associated with any one of a given set of preparation events. We illustrate 
our method by calculating the retrodictive density operator for a two-level 
atom coupled to the electromagnetic field. 



1. Introduction 



Quantum theory is primarily about prediction, that is calculating the probabilities 
associated with any future measurement outcomes given information about the 
preparation, the subsequent evolution of the system and the measurement device to be 
used. A simple example is the probability that an initially excited atom, prepared at 
preparation time t^ , will be found by a suitable detector in its ground state at some 

measurement time t^ . This problem can be solved straightforwardly in the usual 
formalism, which we shall refer to in this paper as the predictive formalism. In this 
formalism, we speak of preparing the system in some state that depends only on the 
outcome of the preparation event. This state is most generally represented by a density 
operator that subsequently changes with time. Associated with the measuring device 
is, in general terms, a probability operator measure (POM) (Helstrom 1976). 
Projection of the density matrix at the time of measurement onto an element of the 
POM determines the probability of the outcome associated with that element. Ideal 
von Neumann measurements of an observable have possible measurement outcomes 
given by its eigenvalues and the associated POM element is simply the projector onto 
the corresponding eigenvector. In this formalism the state represents our knowledge 

of the preparation event. 

Sometimes, for example in communication situations, we only have 
knowledge of the outcome of a measurement event and wish to deduce from this 



information about the preparation event. It is possible to use Bayes' theorem (Box and 
Tiao 1973) together with the predictive formalism for this purpose, but it is more 
direct to use the less common retrodictive formalism of quantum mechanics 
(Aharonov et al 1964, Penfield 1966, Belinfante 1975, Pegg and Bamett 1999, 
Barnett et al 2000a,b). In this formalism a state is assigned solely on the basis of the 

outcome of the measurement. We will refer to such a state as a retrodictive state. 

Thus instead of speaking of preparing the system in a state, we speak of measuring the 
system to be in a state. The retrodictive state transforms as we go back in time to the 
preparation event, at which time projection onto an operator associated with the 
preparation device determines the probability that the system was prepared in the state 
associated with that operator. 

In previous work we have used the retrodictive formalism to study some 
problems in quantum optics including photon antibunching and EPR correlated 
photons, optical state preparation and quantum communications (Pegg and Barnett 
1999, Bamett and Pegg 1999, Bamett et al 2000a,b). In this paper we shall exploit the 
formalism to study some preparation probabilities in a two-level atom interacting with 
an electromagnetic field. 



2. Retrodiction 



2.1. Closed systems 



We refer to a quantum mechanical system as closed if the system is prepared, evolves 



unitarily without interaction with another system and is then measured. A simple 
example involves a spin-half particle which can be prepared in an "up" or "down" 
state with a Stem-Gerlach apparatus. Another such apparatus can give a measurement 
result that the atom is in either of the states |±) corresponding, respectively, to the 
angles 9 and 71 + 6 to the up axis with no change of state between preparation and 
measurement. In one experiment we may know the outcome of the preparation event 
was the up state. We can then use the usual predictive formalism of quantum 
mechanics to calculate the probability that the outcome of the measurement 

corresponds to the state 1+). The result is |(up|-l-)| , that is, cos.^9/2. 

In a second experiment we may be told that the outcome of the measurement 
corresponded to the state 1+) but have no knowledge of the preparation event other 
than that the up and down states had equal a priori probabilities of being prepared. 
We wish to calculate, on the basis of this knowledge, the a posteriori probability that 
the atom was prepared in the down state. For this we would use the retrodictive 
formalism and assign the system the retrodictive state 1+) just prior to the 
measurement. Knowing that the state does not change between preparation and 
measurement we would assign the retrodictive state 1+) to the atom just after the 
preparation event. The probability that the atom was prepared in the down state is 
then K+jdown)! , that is sin^ 0/2. It is customary to speak of the event associated 

with the projection as the collapse of the wavefunction. In the predictive formalism, 
this takes place at the time of the measurement but in the retrodictive formalism it 
coincides with the preparation event. 



We can treat the above simple example more formally in terms of a 
probability operator measure. The essential features of a POM are reviewed in 
Appendix A. For the case above, the measurement POM has two elements, |+)(+| 
and |-)(-| each associated with a possible measurement outcome and which sum to 
the unit operator. The retrodictive density operator associated with a particular 
outcome is just the corresponding POM element, normalised if necessary to have unit 
trace (Barnett et al 2000a). The retrodictive density operator associated with the 
outcome corresponding to detection of state |+) at the time of measurement is 
therefore simply |+)(+|. As there is no evolution in this case between measurement 
and preparation, this is also the retrodictive density operator Pg^'(tp) at the 
preparation time t^ . The source is unbiased because the up and down states have 
equal probability of being produced. The preparation POM elements (see appendix A 
and Barnett et al 2000a) S^ are thus |up)(up| and |down)(down| . The probability 
P(up I +) that the atom was prepared in the up state if detected in the state |+) is 
Tr[pr(^p)S„pj 

Suppose we change the example slightly so that the preparation apparatus 
includes some type of filter, say, which reduces the a priori probability /'(down) for 
the down state to below one-half. Then instead of describing the (now biased) 
preparation device by POM elements S^we use instead the operators A^. These can 

be summed to form the a priori preparation density operator: 



A = Ea, . (2.1) 



p 



In general the terms A^, can be associated with mixed states, but in our example here 
there are two terms, A^p = /'(up)|up)(up| and A^^^^ = P(down)|down)(down|. In the 

absence of any further information as to what state was actually prepared, the best 
description we could give of the prepared state of the atom in the predictive formalism 
would be to assign it the predictive density operator A^p + A^^^^ based on our a priori 
knowledge. If we do know which state was actually prepared, say the up state, then 
we would assign the predictive density operator A^^/TriA\. If, on the other hand, 

we do not know which state was actually prepared but know that the atom is detected 
in the state |+) we can calculate retrodictively the probability that the atom was 
prepared in the up state by projecting the retrodictive density operator corresponding 
to this measurement outcome onto the appropriate term A^ and obtain, after 

normalisation. 



P(^P\+)= ^K\J 'K (2.2) 

£Tr[pr(^p)A,] 



A corresponding expression holds for the probability that the atom was prepared in 
the down state if detected in the state |+) . 

We can summarise the predictive and retrodictive quantum formalisms for 
closed systems as follows. In the predictive formalism A^, associated with outcome 

p of the preparation device, becomes upon normalisation the predictive density 
operator. This operator evolves forward in time until it is projected onto n,„ to 

determine the probability of measurement outcome m. In the retrodictive formalism 
n^, associated with outcome m of the measuring device, becomes upon 
normalisation the retrodictive density operator. This operator evolves backwards in 
time until it is projected onto A to determine the probability of preparation outcome 



2.2. Open systems 



An open system is one in which the subsystem S under study interacts with another 
subsystem. Often this second subsystem is larger and is referred to as the environment 
E. Normally we have at least some statistical information about the initial state of the 
environment and can assign it a density matrix. In the predictive formalism we can 
calculate the probability of obtaining a particular result of a measurement by finding 
the unitarily-evolved predictive density -operator for the combined system at the time 



of measurement and projecting this onto the appropriate combined POM element. In 
the usual case, in which we do not make any measurement of the environment, we can 
write the combined measurement POM element as proportional to the product of the 
POM element n,„ associated with the outcome m for S and the unit operator Ig on the 

space of the environment states. An alternative predictive approach is to represent the 
subsystem 5 by a reduced density operator at the time of preparation, which is the 
trace of the density operator of the complete system over the states of the 
environment. The subsequent non-unitary evolution of the reduced density matrix is 
determined by a master equation. This equation is constructed so as to produce the 
correct measurement probability when the evolved reduced density operator at the 
time of the measurement is projected onto 11^ . Thus, although the complete system 
usually evolves into an entangled state we can, in the usual predictive formalism, 
represent the state of S in the limited sense that from it we can calculate the 
probability for S to be detected in any particular state, conditioned on no information 
being obtained about the state of the environment. 

In the retrodictive formalism of quantum mechanics, the retrodictive density 
operator of the system being measured is, immediately prior to the measurement, 
simply the POM element associated with the measurement outcome with a suitable 
normalisation factor (Barnett et al 2000a). Thus if the measurement gives an outcome 
m for S but no information about the environment, we can find the retrodictive density 



operator at an earlier time t as 



^».(,^)_ ii%Kv 



Tr, 



ES 



IeH. 



(2.3) 



where Uh the forward time shift operator from the time t^ to time t^ , the time of the 

measurement. 

The terms A^^ of the preparation device density operator for the system 

comprising S and E can be written as A^ s^pE- If the preparation device is set up 
such that in the predictive formalism the environment always has the same initial state 
p^^^itp), then A^ = A^^ sPe"'''(^p) • The probability that the system was prepared at 
time t in the state corresponding to A 5 is thus given by 



Tr, 



ES 



^...^ ,, 1 Tr,4^-1,fl„^A,,,pr(g- 



Trs[prA^,, 



lTr,[prA^,, 



(2.4) 



where 



■Iretr 

Fs ~ 



Tr, 



pr\t^)u'hj 



Tr,,[pr'%)U'U^U] 



(2.5) 



The reduced retrodictive density operator Pg"" is of course not the same as the 
reduced predictive density operator. Nevertheless it is "reduced" in the sense that it 
acts only on the state space of S and not on the environment states. As with the 
reduced predictive density operator, this operator can be used to calculate only 
probabilities conditioned on the measurement not providing any information about the 
environment. 

The main aim of this paper is to calculate Pg"" , which we shall refer to simply 

as the retrodictive density operator, for some open atomic systems. We can obtain 
the retrodictive density matrix elements in terms of a complete orthonormal basis set 
of states {|/)} of S. We have from (2.5) 



{i\pr\k) = -Tv,[{i\pr(t^)Wn^U\k) 



— Trps 



flAmipr'iOU'] (2.6) 



where we have used the cyclic property of the trace and A'^ is the denominator in (2.5). 
By introducing a general time-independent operator in the form 



As = Y.AMk\, (2.7) 



10 



we can rewrite the retrodictive density matrix elements (2.6) as 



{i\pr\k) = ^^Tr,,[hMpr\)U'] ■ (2.8) 



If A^ has the form of an initial density matrix then we can find 

Trgs[n^t/AsPg"' (?p)f/ ] from the appropriate master equation. We can most easily 

find the normalisation constant A'^ by noting that 



K^-^rio = ^Y,j^^^AKuKpT\Km = i- (2.9) 



Expression (2.8) is useful for obtaining the retrodictive density matrix elements 
directly without first calculating preparation probabilities from (2.4). The retrodictive 
density matrix can then be used to calculate the preparation probability for any state 

by projection onto the associated preparation POM element S^ or preparation density 

operator term A^, as appropriate. 

An alternative approach, if we wish to calculate the preparation probability 
associated with a particular preparation POM element S^ without explicitly finding 

the retrodictive density operator, is as follows. This probability is given by 
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Tr^prS,) = ^Tr,,[pr\t^)U%pZ^]. (2.10) 



We can rewrite this using the cyclic property of the trace as 



(prsJ = i-TrJn,„S/^j], (2.11) 



Trc, 



where 



Z^(tJ = Tv,UZ^(t^)pr(t,)U' . (2.12) 



If S^ has the form of an initial density matrix we see from (2.12) that S^,(?^) has the 

form of a reduced density matrix in the predictive formalism obtainable from the 
appropriate master equation for the initial state of the environment (Barnett and 
Radmore 1997). It is straightforward to generalise this approach for a preparation 
density operator term A^ in place of S^ . 

In general, finding the retrodictive density operator is more useful in that it 
allows the subsequent calculation of the preparation probability of any state and this is 
the approach we use for most of this paper. We do, however, illustrate the alternative 
approach of the above paragraph by finding directly the preparation state probabilities 
for a spontaneously emitting atom. For completeness, we also outline in Appendix B 

12 



another method devised specifically for calculating the retrodictive density operator 
for a two-level atom. 



3. Two-level atomic retrodiction 



The retrodictive formalism for open systems allows us to calculate, from the outcome 
of a measurement, the probability associated with an earlier preparation event. In the 
preceding section we have presented formulae for obtaining directly either the 
probabilities for particular state preparations, or the retrodictive density matrix 
elements in any basis, which can then be used to calculate preparation probabilities. 
Here we shall illustrate our methods by using some examples from atomic physics, all 
based on the two-level atom coupled to the electromagnetic field. We work 
throughout in an interaction picture in which the energies of the two atomic levels do 
not appear in the Hamiltonian. 



3.1. Spontaneously emitting atom 



Consider first the simplest of open systems, the damped undriven two-level atom with 
excited and ground states |e) and |g). We assume that a measurement has been made 
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at time t^ and that we know the state corresponding to the measurement result. We 
wish to retrodict a preparation event that has taken place at an earlier time t^=t^-x. 

Equation (2. 1 1) gives the probability that the atom was prepared in a state 
corresponding to a particular preparation POM element. We simply need the 
appropriate measurement POM element and the precise form of 'Ep{t^). To be 
specific, suppose we find the atom to be in the ground state, so that fl^ = |g)(g| is our 
measurement POM element. We seek, respectively, the probabilities that the atom 
was prepared in the excited state and the ground state, given that one of these two 
states was prepared at time t^ by an unbiased apparatus so that the a priori 

probabilities for these two preparation events are 1/2. For a measurement made 
immediately after preparation ( T = ) the ground state preparation probability should 
be unity and the excited state preparation probability should vanish. In the other 
extreme, if there is a very long time between preparation and measurement, detecting 
the atom in the ground state should give us no information at all about the preparation 
state as all states decay eventually to the ground state. Thus the ground and excited 
state preparation probabilities given that we observe the atom in its ground state are 
1/2 and the measurement has given us no information about the preparation of the 
atom. We now find the probability that the atom was prepared in the excited state at a 
general time T prior to the measurement for which the corresponding preparation 
POM element is E^^ = |e)(e|. This element has the form of a density operator so we 
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can find S (?^) from the usual predictive master equation (Bamett and Radmore 



1997) obtained by taking the initial state of the environment pl^^^itp) as the vacuum 



state. The solution to this master equation is 



Sp(C) = |e>(e|exp(-2rT) + |g)(g|[l-exp(-2rT)], (3.1) 



where T is one half of the Einstein A-coefficient for spontaneous emission. The 
excited state preparation probability conditioned on us detecting the atom in the 
ground state is then 

pP^'^-'Ce I g) = ^Tr4|g)(g|[|e)(e|exp(-2rT) + |g)(g|(l - exp(-2rT))]} 

= ^[l-exp(-2rT)] . (3.2) 



Similarly, the probability that the atom was prepared in the ground state can be 
calculated using S^^ = |g)(g|. This does not evolve with time as an undriven atom 

initially in the ground state must remain there. The ground state preparation 
probability thus becomes 



^""(glg) = ^Tr,(|g)(g|g)(g|) = i^ . (3.3) 
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The probabilities (3.2) and (3.3) must sum to unity and hence 



A^ = 2-exp(-2rT), (3.4) 

which gives 



jprep. ,„^_ [1 - exp(-2rT)] 



pp-p(e I g) = ± fi^ ^ (3.5) 

[2-exp(-2rT)] 



P''^'=''(glg) = ^ . (3.6) 

[2-exp(-2rT)] 



Clearly for T = , P'^'^e I g) = and P'^'^g I g) = 1 as anticipated above. For 
T ^ oo , P''"'''(e I g) = P''"'''(g I g) = 1 / 2 , so we have no preparation information. A 
similar calculation gives the excited and ground state preparation probabilities 
conditioned on us detecting the atom in the excited state. The results, P''"'''(e I e) = 1 
and P''"'''(g I e) = for any time interval between preparation and measurement, are 
intuitively obvious. If we detect the atom in the excited state then it cannot have 
decayed and it must, therefore, have been prepared in the excited state. 

It is also possible to prepare and measure the atom in superpositions of the 
excited and ground states. Any measured superposition containing a non-zero ground 
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state component gives no information about the prepared state if a long time has 
elapsed between preparation and measurement. 

The above method for retrodicting probabilities required us to specify the 
preparation POM element in advance. A more versatile method is first to calculate 
the retrodictive density matrix associated with a measurement outcome. This density 
matrix can then be projected onto any preparation POM element. The retrodictive 
density matrix elements are readily obtained from equations (2.8) and (2.9). We write 



^s=Ae|e)(e| + A^Jg)(g| + Aje)(g| + A^Jg)(e|, (3.7) 



and let this have the properties of a density operator. The appropriate master equation 
can then be solved to give 



= \e){e\A^^ exp(-2rT) + \g){g\{A^, + 4J1 - exp(-2rT)]} 

+|e)(g|4^exp(-rT) + |g)(e|A^,exp(-rT) . (3.8) 



The retrodictive density matrix elements can now be found from (2.8). These 
elements are 
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(e|pr|e) = ^^Tr,[A,ajn„] 

= ^ {exp(-2rT)(e |n„ I e) + [1 - exp(-2rT)](g |n„ | g)} , (3.9) 



(g|pr|g) = ^(g|njg), (3.10) 



(g|pr|e> = ^(g|nje)exp(-rT) (3.11) 



and 



(e|pr|g) = ^(e|nJg)exp(-rT) . (3.12) 



The nonnalisation constant is such that the sum of the diagonal elements is unity. 
Setting the sum of (3.9) and (3.10) to unity gives 



iV = exp(-2rT)(e|nje) + [2-exp(-2rT)](g|nJg). (3.13) 



We note that for T = 0, the retrodictive density matrix is simply Il^/Tr(n^) . This 
simply states that the retrodictive state immediately before the measurement is 
proportional to the POM element associated with the measurement outcome (Barnett 
et al. 2000a). 
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As an illustration, we use this retrodictive density matrix to calculate the 
probability P'"^''''(e I g) that the excited state was prepared by an unbiased apparatus 
which is capable of preparing the atom in the excited or ground states, if the atom is 
detected to be in its ground state. Projecting the retrodictive density matrix onto the 
preparation POM element |e)(e| gives for this probability the expression (3.9) with 
n^ = |g)(g|. This gives the result (3.5) obtained previously by direct calculation. 

The retrodictive density matrix gives us more than just the ground and excited 
state probabilities. As an example, suppose that our measurement of the atom gave a 
result corresponding to the superposition state \6) = (cos(0/2)|g) + sin(0/2)|e)) . 

Furthermore, let us suppose that we know that the atom was prepared in one of the 
non-orthogonal states |e) and 1+) = (|e) + |g))/A/2 and that the prior probabilities for 

these two states are p and I- p respectively. How does our measurement result 
change these probabilities? To answer this question we first calculate the retrodictive 
density matrix associated with our measurement result for which the POM element is 



n. = |W| = (cos(0/2)|g) + sin(0/2)|e))(cos(0/2)(g| + sin(0/2)(e|). (3.14) 



A little algebra then shows the retrodictive density matrix elements to be 



19 



(e|prie) 



- {l + cos 0[l - 2 exp(-2rT)]} 
l + cos0[l-exp(-2rT)] ' 



(3.15) 



1 



(g|prig)= ^ 



(1 + COS0) 



l + cos0[l-exp(-2rT)] 



(3.16) 



and 



(g|prie>= 



sin0exp(-rT) 



l + cos0[l-exp(-2rT)] 



= (e|pr|g) 



(3.17) 



The two possible preparation events are associated with the operators 



Ae=p|e>(e| 



(3.18) 



and 



A. = (l-p)|+)(+l = ^(l-p)(|e> + |g))((e| + (g|). 



(3.19) 



The probabilities that the atom was prepared in the state |e) or |+) given that it was 



subsequently observed to be in the state \0) are 



Tr. 



pP"''(e I 6) 



pTK 



Tr. 



pr(K+K) 



p{l + cos 0[l - 2 exp(-2rT)]} 



1 + cos 0[l - exp(-2rT)] + (l- p)s'm6 exp(-rT) - p cos 6 exp(-2rT) 



(3.20) 
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and 



pp-p(+ 10) = 



Tr< 



prK] 



Trc 



pr(Ae+A.)] 



(1 - p){l + COS 0[l - exp(-2rT)] + sin exp(-2rT)} 
1 + cos 0[l - exp(-2rT)] + (l- p)smO exp(-rT) - p cos exp(-2rT) 



(3.21) 



For T = , these probabilities are 



pP'^Ce I 6>) 



p(l-cos0) 



p(l - COS 0) + (l- p)(l + sin 0) 



(3.22) 



pP'^P(+ I 0) = 



(l-;?)(l + sin0) 



p(l - cos 0) + (1 - p)(l + sin 0) 



(3.23) 



These values agree with those obtained using Bayes' theorem directly (Bamett et al 



2000a), 



/'P"''(e I 0) = 



P(e,6>) ^ p|(e|6>)| 

P(0) p\{e\0)f+il-p)\{+\0)f 



(3.24) 



,prep,^m^_^(+'0)_ (1-P)|(+|0) 



pprep(^ I 6>) 



Pi0) p\{e\0)\+(l-p)\{+\0)\ 



2 ■ 



(3.25) 
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For T ^ oo , the probabilities reduce to pP'^^Ce 10) = ;? and pP''=P(+ \e) = l-p, that is, 



the prior probabilities. This tells us that, after a long time, the measurement gives us 



no additional information about the initially prepared state. 



3.2. Incoherently excited atom 



We now consider the case of our two-level atom coupled to a thermal field of mean 
occupation number n . This modifies the master equation governing the evolution of 



the density matrix (Barnett and Radmore 1997) and in place of (3.8) we obtain 



MtJ = 



|e)(e| 4, 



— ^^ + exp(-2r(2n + 1)t) ^^-^ + A^^ -^^— [l - exp(-2r(2n + 1)t)1 



+|g)(g|4 



— — + -^exp(-2r(2n + 1)t) 
,2n+l 2n+l 



+ 4e ^^[1 - exp(-2r(2n + 1)t)]J 



+1 e>(g |4„ exp(-2r(2n + 1)t) + 1 g)(e U, exp(-2r(2n + 1)t) . 



(3.26) 



The retrodictive density matrix elements can again be found by differentiation using 



(2.8). This procedure gives 
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(e|pr|e> = ^ (einje; 



" +-^^exp(-2r(2n+l)T) 



2n+l 2n + l 



n +1 



+(g|n Jg>^3:^[l - exp(-2r(2n + 1)t)] 



(3.27) 



(g|prig)=- (g|njg; 



|:^ + -^exp(-2r(2n+l)T) 
2n+l 2n+l 



+(e|n Je)-^[1 - exp(-2r(2n + 1)t)] , (3.28) 

2n +1 J 



g|pr|e) = -(g|nje)exp(-r(2n+l)T) 



A^ 



(3.29) 



and 



e|pr|g) = -(e|nJg)exp(-r(2n+l)T). 



A^ 



(3.30) 



The normalisation constant is 



A^ = (e|n^|e: 



2n \ 

— — + — — exp(-2r(2n + l)T) 
2n +1 2n +1 



+{s\A„\g] 



2n+2 1 



2n + 1 2n + 1 



exp(-2r(2n + 1)t) 



(3.31) 



For T ^ oo , it is clear that the retrodictive density operator tends to one half of the 



identity operator, corresponding to no information about the initially prepared state. 



This happens for all measured states because the thermal radiation can cause 
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absorption so that it is possible to prepare the atom in its ground state and to find it 
later in its excited state. Figure 1 shows the predictive and retrodictive diagonal 
density matrix elements for an atom respectively prepared or measured to be in its 
excited state. The predictive density matrix elements rapidly decay to their steady 
state form for which the ratio of the excited and ground state probabilities is 
n/(l + n) . The retrodictive diagonal density matrix elements tend to the value 1/2 for 
long time intervals T . 



3.3. Coherently driven atom 



In our final example the atom is resonantly driven by an intense laser field. It is well 
known that such a system exhibits Rabi oscillations at a frequency 
Q. = (V^ - r^ /4)"^ where Vis proportional to the laser field strength and T is the 
atomic spontaneous decay rate. The (interaction picture) Hamiltonian describing the 
coupling between the atom and the laser field is 



H = ^{a,+a_), (3.32) 



where d^ = |e)(g| and a_ = |g)(e|. The reduced, predictive density operator for the 
atom satisfies the master equation 
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dt 



pr = -i 



V 



{a^ + a_),p. 



pred 
S 



+ r[2a_pr'a^ - a^ajr' - pr'^,a_]. (3.33) 



The solution for the predictive reduced density operator at the measurement time may 



be written in the form 



^ pred I 



Ir 



pr (O = 2^ + "(^)^i + ^(^)^2 + w(T)(y, 



(3.34) 



where 1 = |e)(e| + |g)(g| is the unit operator on the space of the atom. The Pauli 
operators di = |e)(g| + |g)(e|, dj = /(|g)(e| - |e)(g|) and dj = |e)(e| - |g)(g| represent 



respectively the real and imaginary parts of the dipole moment and the atomic 



inversion. The functions m(t), v(T)and w(t) are the expectation values of dp dj and 



a^ respectively at time T after preparation and have the values (Barnett and Radmore 



1997) 



M(T) = M(0)exp(-rT) 



(3.35) 



v(T) = —^ — ^(irv + exp(-3rT/ 2){ [(y' + 2r')v(0) - 2ry]cos(^T) 



+ 



a"'[-3rV + (y' + 2r')(-yM'(0) + rv(0)/2)]sinaT}) 



(3.36) 



W(T): 



y" + 2r 



- — ^(-2r' + exp(-3rT/2){[2r' + (y' + 2r')w(0)]cosaT 



2 , oT-2 
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+a-'[-2r(v^ + r' /2) + (v^ + 2r')(yv(0) - rw(0)/2)]sinaT}). (3.37) 



The excited and ground state matrix elements can be found from the master equation 
solution (3.34) using 



e|Ps (?p)|e> = 4e =[l + >v(0)]/2 (3.38) 



g|Ps (?p)|g> = A,g =[l-w(0)]/2 (3.39) 



e|Pr'(^p)|g> = 4g = ["(0) - iv(0)]/2 (3.40) 



g|Ps (?p)|e> = 4e = [«(0) + M0)]/2. (3.41) 



These can be substituted into (3.34)-(3.37) which allows us to calculate the 
retrodictive density matrix elements according to (2.8). The calculations are 
somewhat lengthy to reproduce in detail, so we shall concentrate on the results for a 
few specific measurement POM elements. 
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We first study the case where the atom is detected to be in its excited state, so 
the measurement POM element is ft^ = |e)(e|. The retrodictive matrix elements are 



then 



(e|pr|e) = 



2N 



^ I dw{t) 

V ^4e J 



1 



2NiV' + 2Y') 



-{V + exp(-3rT/2)[(y' + 4r')cos^T 



20. 



(5y'+4r')sinaT]}, 



(3.42) 



(g|Ps |g) = 



V' 



2N{V^+2T^) 



( 3r 

l-exp(-3rT/2) cos(Qt)H sin(QT) 

V 2Q. 



(3.43) 



-iV 



(e|pr |g) = -(g|pr |e> = -^exp(-3rT/2)sin(^T), 



2NQ. 



(3.44) 



where the normalisation constant is 1 + 



dw{t) dw{t) 
dA dA. 



. The evolution of the 



gg J 



predictive and retrodictive matrix elements for this system is shown in figure 2. In the 
retrodictive case the diagonal matrix elements show Rabi oscillations which decay to 
the "no information" level of 1/2. The off-diagonal elements decay to zero as 
expected. This corresponds to a loss of coherence information about the system a 
long time in the past. Similar behaviour is found for n^ = |g)(g|. 

Suppose instead that the atom is measured in the superposition ^=(|e) + z|g)) , 
which gives n„ = — (1 + Cj). The retrodictive matrix elements are 



(e|pr|e) = -^ 1 + 



2N \ V + 2Y' 



2YV 
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-exp(-3rT/2)[2rVcos(nT) + — (y'+5r')sin(aT) 



(3.45) 



(g|prig)=^ 1+^^ 



irv 



-exp(-3rT / 2)( 2r Vcos(aT) + — (r' - y' ) sm{Q.r) 



(3.46) 



(e|Ps |g) = -(g|Ps |e) = 



-/exp(-3rT/2) 



2N 



cos(Qt) h sin(QT) 

2a 



(3.47) 



with A^ = 1 + 



(9v(t) (9v(t) 



. These are shown in figure 3, along with the 



corresponding predictive evolution. At the measurement time the retrodictive diagonal 
matrix elements take their long-time values of 1/2. The off-diagonal matrix elements 
are pure imaginary, and so they undergo Rabi oscillations. These elements couple to 
the diagonal matrix elements so that we also see oscillations in the diagonal density 



matrix elements. 



Suppose now that the atom is measured to be in the superposition 
-^(|e)-l-|g)), so that the measurement POM element, ft,,, = —(l -I- di), contains only 



the real part of the dipole matrix element. In this case the retrodictive matrix elements 



are much simpler. 



(e|pr|e) = (g|pr|g) = 



(e|pr|g) = (g|Pr|e> = ^exp(-rT). 



(3.48) 
(3.49) 
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Here the diagonal matrix elements are time-independent. The reason is that the off- 
diagonal matrix elements are purely real and so do not lead to Rabi oscillations, that is 
m(t) simply decays. 

The predictive state of a driven two-level atom evolves forward in time 
towards a steady state with an unchanging population inversion and imaginary dipole 
matrix element. The expectation values of the Pauli operators a^, d^ and O-^ for this 



state are 



u{x)-- 


= 0, 




V(T) = 




lYV 


~ V 


^ + 2r^' 


^„(^\ 




-2r' 



(3.50) 
(3.51) 



(3.52) 

v + 2r' 



If we were to perform a measurement showing the atom to be in this state then we 
could calculate a corresponding retrodictive density operator. The diagonal and off- 
diagonal matrix elements of this density operator are shown in figures 4a and 4b 
respectively. The off-diagonal matrix elements have a non-zero imaginary part so the 
retrodictive density matrix exhibits damped Rabi oscillations. The long time interval 
limit of the retrodictive density operator is, once again, one half of the identity 
operator, corresponding to providing no information about the initially prepared state. 
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4. Conclusions 



The retrodictive formalism provides the means to evaluate preparation 
probabilities given the results of some later measurement event. Our approach relies 
on the introduction of a retrodictive state. This is a state assigned to the system prior 
to the measurement, based on the measurement result. Immediately before the 
measurement, this retrodictive state is simply the normalised POM element 
corresponding to the measurement result (Barnett et al 2000a). Put more simply, the 
retrodictive density operator is the state corresponding to the measurement result. If 
the system is closed, that is isolated from its environment, then the change in time of 
the retrodictive density operator is unitary. If, however, the system is coupled to an 
unmonitored environment then we will lose information about the system through 
coupling to the environment. This will mean in most cases that the retrodictive density 
operator for long times prior to the measurement will tend to become proportional to 
the identity operator. This is the zero information state and reflects the fact that our 
measurement tells us nothing about the way the system was prepared long before the 
measurement. 

We have applied the formalism to retrodict from measured states of a two- 
level atom interacting with the electromagnetic field. For an undriven atom, which 
can only undergo spontaneous emission, we find that the diagonal elements tend to 
the value 1/2 and the off-diagonal ones to zero. This reflects the fact that we can 
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obtain from the measurement no information about the system a long time before the 
measurement. The only exception is if the atom is measured in its excited state. In this 
case, the retrodictive density operator must be the excited state projector |e)(e| for all 
times prior to the measurement. This is because there is no light to absorb and so the 
atom could not have made a transition to the excited state from the ground state. With 
this exception, the decay to "no information" is a general feature of retrodiction for 
damped two-level atoms, both driven and undriven. 

We have also performed retrodictive calculations for the incoherently and 
coherently driven atom. The evolution of the retrodictive density matrices of both 
systems is different in character from that of their predictive counterparts. In the 
forward time direction an atom driven by an incoherent field relaxes to a steady state 
predictive density operator with the ratio of the excited and ground state probabilities 
of n l{n + 1). The corresponding ratio for retrodictive evolution is unity. For the 
coherently driven atom we find retrodictive Rabi oscillations which decay. Again the 
predictive and retrodictive evolutions decay to different forms for long time intervals, 
with the retrodictive state tending to the zero-information state with density operator 

i/2. 
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Appendix A. Probability operator measures 



A comprehensive discussion of a probability operator measure (POM) is given by 
Helstrom (1976). Here we outline only the essential features. 

In the predictive formalism we can use a POM as a mathematical 
representation of a measuring device. The POM is a set of elements n„, , each of 
which is a positive semi-definite operator associated with a different possible outcome 
of a measurement. n„, is defined such that the probability of the outcome m from a 
measurement of a system with a predictive density matrix p''""' at the time of 
measurement is given by Tr[p''""'n^] . The total probability for all possible outcomes 
must be unity. Thus 



J^Tr[p-^nj=Tr[p''-'^£nj 

m m 

= 1. (A.l) 



For this to be true for all p""^*"^, normalisation of the trace of p''""' to unity requires 
^ n^ to equal the unit operator. 

m 

A measurement that yields no information may be considered as a 
measurement with only one possible outcome, for example a zero meter reading, for 
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all p'"'"'. Normalisation of the POM as above requires the corresponding single 
element to be the unit operator. 

We can also use a POM with elements E^ as a mathematical description of an 

unbiased preparation device. In this case Tr[p^''"Sp] is the a posteriori probability of 

a preparation outcome p based on knowledge of the measurement outcome. This 
probability is equal to that obtainable from the predictive formalism combined with 
Bayes' theorem (Barnett et al 2000a). If the preparation device is biased we work 
instead with components A^ of the a priori density matrix. 
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Appendix B. Alternative derivation of retrodictive density matrix 



Here we outline an alternative method for finding the retrodictive density matrix for a 
two-level atom coupled to the environment. We first expand the retrodictive density 
operator for the atom in terms Pauli operators Oj with j = 1,2,3 (Barnett and 
Radmore 1997). The general form of this expansion is 



•iretr _ ^ 

Ps -2 



r 



V j J 



(B.l) 



where 1 is the unit operator on the space of the atom. Multiplying by (1 + d^) and 
tracing over the atom states we obtain 



Tx,[{\ + a,)pT] = l + u, . (B.2) 



Thus from (2.5) we have 






u, = ^TY^,V{\ + a,)p^^{t^)U'YiJJ]-\ 



= ^Tr,,[i7(i + (7jpr(fp)^'nj-l 



= ^Tr,[(i + <7jajnj-l (B.3) 



where, because (1 + (7j)/2 has the form of a density matrix, we can find 
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(i + a,)(tj = Tr^[U(l + a,)pr\)U'] 



(B.4) 



from the appropriate master equation. Indeed the values of (1 + a,^)(t^) can be 



obtained immediately from the general expressions given by Barnett and Radmore 
(1997), that is, equations (3.35)-(3.37) with values of m(0), v(0) and >v(0) chosen as 



or 1 as appropriate. Inserting these into (B.l) gives the retrodictive density matrix 



in the form 



pr = 



l + l{^Tr,[(i + (Tp(.jnj-l}d^. 



(B.5) 



Finally we note that the expression for A^, which is given by the denominator of (2.5) 



can be written as 



N = Tr^[Upr'(t,)®lu'flJ 



= Tr,[l(rjnj 



(B.6) 



where, because 1/2 has the form of a density operator, we can find Kt^) = 



Tr, 



^ A, 1 "^ ^ -I- 1 

Up^" (t ) ® It/ from the same master equation as used above. 
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Figure Captions 

Fig. 1 . Incoherently excited atom with a mean thermal occupation number of n = 1 . 
(a) Predictive diagonal density matrix elements plotted against time for an atom 
initially prepared in the excited state, (b) Retrodictive diagonal density matrix 
elements for an atom measured to be in the excited state. The solid and dashed curves 
are the excited and ground state matrix elements respectively. The time T, measured in 
units of the atomic state lifetime, is the interval between preparation and 
measurement. Note that in all figures, the preparation time is to the left of the 
measurement time. This means that retrodictive density matrix elements are plotted 
with T increasing to the left. 



Fig. 2. Coherently driven atom with V = 4r . (a) Predictive diagonal density matrix 
elements for an initially excited atom, and (b) retrodictive diagonal density matrix 
elements for an atom measured to be in the excited state. The solid and dashed curves 
are the excited and ground state matrix elements respectively, (c) Imaginary part of 
the predictive and (d) retrodictive off-diagonal matrix elements for the same systems 
as (a) and (b). The solid and dashed curves are the imaginary parts of (elPs^^jg) and 
(g|Pr''l^) respectively. The real parts are zero. 
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Fig. 3. Coherently driven atom with V = 4T . (a) Predictive diagonal density matrix 
elements for an atom initially in the superposition (|e) + i\g)y-\f2 , and (b) retrodictive 

diagonal density matrix elements for an atom measured to be in the same 
superposition, (c) Imaginary part of the predictive and (d) retrodictive off-diagonal 
matrix elements for the same systems as (a) and (b). The solid and dashed curves 
represent the same quantities as in figure 2. 



Fig. 4. Coherently driven atom with V = 4r . (a) Retrodictive diagonal density matrix 
elements for an atom measured to be in the predictive steady state, (b) Retrodictive 
off-diagonal matrix elements for the same system. The solid and dashed curves 
represent the same quantities as in figures 2b and 2d. 



39 








0.2 



0.4 ,^0.6 

x/r 



0.8 



1 




1 



0.8 



0.6 ,^0.4 

x/r 



0.2 







0.4" 



-p 

Q) 

e 

Q) 



0.2" 



H 



/ N 



o 

"^-0.2 
I 

M-l 
M-l 
O 

-0.4 + 




/ \ "0.4 



5 




0.8t 




tO.8 



' \ |0 . 7 

o 

H- 

■0.6^ 




0.4" 



-p 

Q) 

e 

Q) 



0.2" 



H 



0' 



o 

"^-0.2 
I 

M-l 
M-l 
O 

-0.4 + 



^ X 



1 



2 3 

T/r 



5 



5 



I ■ ■!■ 



3 



x/r 





\"0 . 6 



o 

H- 

O 

0) 



--0.5 



I— ' 
CD 

3 

(D 

rt 



"0.4 



5 






3 




